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We consider propagation of gravitational radiation in a magnetized multicomponent plasma. It 
is shown that large density perturbations can be generated, even for small deviations from flat 
space, provided the cyclotron frequency is much larger than the plasma frequency. Furthermore, 
the induced density gradients can generate frequency conversion of electromagnetic radiation, which 
may give rise to indirect observational effect of the gravitational waves. 
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I. INTRODUCTION 

Recently there has been an increased interest in gravi- 
tational waves, mainly due to possibility of direct detec- 
tion by LIGO (Laser Interferometer Gravitational- Wave 
Observatory) Naturally the effects of gravitational 
waves on earth are very small - which is illustrated by 
the large dimensions required for detection. Closer to 
the source the influence of the gravitational waves may 
be larger, but generally it is nontrivial to predict the 
possible influence of the emitted radiation - in particular 
the coupling to the electromagnetic fleld complicates the 
description. For a discussion of the interaction between 
electromagnetic flelds and gravitational radiation in an 
astrophysical context, see for example Refs. and 
references therein. 

In the present paper we will study the propagation 
of gravitational perturbations in a magnetized plasma, 
with the direction of propagation perpendicular to the 
magnetic fleld. It turns out that large density gradients 
driven by the gravitational perturbation can be gener- 
ated, even for small deviations from flat space, provided 
the cyclotron frequency is much larger than the plasma 
frequency. Furthermore, as is well known from labora- 
tory plasmas (see, e.g., |^), moving density gradients can 
increase (or decrease) the frequency of electromagnetic 
wave packets, so called photon acceleration. The den- 
sity gradients in our case are propagating with exactly 
the speed of light, in contrast to the laboratory appli- 
cation H]. In principle this means that a given photon 
may increase its energy by several orders of magnitude, 
independent of its initial energy. Applying our results 
to gravitational radiation generated by binary systems, 
it turns out that the regime of most interest is the in- 
frared regime. In this case frequency conversion by an 
order of magnitude is possible, for a binary system close 
to merging. 



II. PLASMA RESPONSE TO A GRAVITATIONAL 
WAVE PULSE 



A. Basic equations 

The metric of a linearized gravitational wave propa- 
gating in the z-direction can be written as 

ds^ = -dt^ + [1 + h{u)]dx^ + [1 - h{u)]dy^ + dz^ , 

(1) 

where we have assumed linear polarization, and u = z — t. 
For an observer comoving with the time coordinate, the 
natural frame for measurements is given by 

eo = 9t , ei = (l - \h) , 62 = (l + \h) dy , 63 = 9^ . 

(2) 

It can be shown |0] that in such a frame, Maxwell's equa- 
tions can be written 



dE 

'dt 

dB 

'dt 



V E = p/eo , 
V-B = , 

- VxB = -j^-noj , 
+ VxE = -j^ , 



(3a) 
(3b) 

(3c) 
(3d) 



where the effective gravitational current densities are de- 
flned as 



and V = (61,62,63). 

To flrst order in h, the fluid equations become 

Qfi 

— + V-{nv) = 0, 



d_ 

dt 



q 



v-V "fv ^ —{E + vxB 



(4a) 
(4b) 

(5a) 
(5b) 



where 7 = (1 — i;|) v\\ = V3, and n — jn, where 
n is the proper number density. These equations holds 
for each particle species. Note that in general terms pro- 
portional to Vih and appear in the equations [||. 
Throughout this paper, we will assume that vi,V2 ^ 1, 
and thus neglect terms of order vih, V2h. 
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B. Electromagnetic fields driven by a gravitational 
perturbation 

Prom now on we assume d/dt ^ Wc, where Wc = qB/m 
is the cyclotron frequency, for all particle species (since 
the gravitational perturbation is assumed to be the driver 
of all perturbations this scaling thereby holds for d/dt 
acting on all fields). Furthermore, we assume the pres- 
ence of an external magnetic field: Bq = BqCi [where the 
total field is B = [Bq + SB)ei]. The electric field takes 
the form E ~ E±e2. 

Looking for solutions driven by the gravitational per- 
turbation, and thus using d/dt = —d/dz, we first con- 
sider Faraday's law for SB <^ Bq, which gives 



6B = -Ej_ + hBo 



(6) 



Next we note that if the excited fields E±_ and SB grows 
(invalidating SB <ti Bq), the quantity E± + B that ap- 
pears in the effective current still becomes E± + SB = 
hBo, and thereby the above formula holds for arbitrary 
electromagnetic amplitude. Taking the time derivative 
of Ampere's law, using Eq. (H), we obtain 



dt^ 



d^ 

dz^ 



E, 



Mo 



dt 



1^2 Bo , (7) 



where the sum is over particle species, and j± = j2- For 
d/dt = —d/dz, the term (explicitly) involving Ei_ van- 
ishes. The currents are determined by the equation of 
motion, noting that the condition d/dt <^ uJc means that 
the current contribution from different particle species 
cancel to lowest order in an expansion in the operator 



^d/dt. The equation of motion gives 

Ei_ 



Bo + SB 



(8) 



to lowest order. Note that, using Eq. (0), we can now 
approximate the denominator in Eq. (g) by Bq — E±. 
The error this approximation introduces will not have 
any noticeable effects. This is because can only be 
altered significantly by the omitted term if SB w Bq, 
but this regime is inaccessible, since - from Eq. (||) - 
it correspond to superluminal speeds. From the parallel 
component of Eq. (5b) we can calculate the first order 
correction to the induced velocity, which subsequently 
determines the current. We obtain 



1 — W|| d{'jv\\ ) 



q Bo- E± dt 
Furthermore, the continuity equation gives 

noV\\ 



Sn 



1 



(9) 



(10) 



where we have divided the density into a perturbed and 
an unperturbed part, n = no + Sn. 



>From (0) and the relations above we can thus de- 
termine the induced velocity and density in terms of the 
metric perturbation h. The result (for all particle species) 



i-{i-ny 



1 + {1-H) 



Sn=- 



1 



- 1 



(11a) 
(lib) 



where H 



Pii) 



and 



'Mi) 



(9(^i)"'o/£o"^(j))^^^ is the plasma frequency for the un- 
perturbed plasma species i. Thus it is clear that even 
a moderate or small value of the gravitational pertur- 
bation may cause significant density perturbation, pro- 
vided the plasma is strongly magnetized in the sense 
that X)i('-^p(i)/'^c(i)) ^ ^- This is because the fast 
magnetosonic (or compressional Alfven) wave fulfills ap- 
proximately the same dispersion relation as the grav- 
itational wave, with the mismatch being proportional 
to J2ii^p(i)/^c{i)) 101 • '^^^ divergence that occur for 
7i 1 is clearly unphysical, and it's removal will be 
discussed in the next subsection. 

For future considerations it will also be useful to have 
the relation between the relative magnetic field pertur- 
bation and the relative density perturbation. When 
\SB\ ^ |ft-Bo|, which is the case of most interest, the 
last term of Eq. (^ can be neglected and the desired re- 
lation can be derived by combining the resulting formula 
with Eqs. (||) and ([lo|). The simple result is 



Sn 
no 



S_B_ 

B^ 



(12) 



C. Removal of the divergence 



The purpose in this subsection is to explain the rea- 
son for the occurrence of divergence when Ti. approaches 
unity, and to discuss various modifications of the assump- 
tions that lead to a more physical behavior. From Eq. 
(^ we note that for infinitesimal velocity perturbations, 
v± (and thereby j±) depends linearly on v\\, but for large 
parallel velocities, in particular when v\\ ^ 1, v± remains 
finite due to the factor 1 — W|| . From Eq. (Q) it is thus clear 
that we cannot have a stationary solution where E \ de - 
pend only on z — t for large enough h, and from Eq. (11a) 



we see that this limit for the gravitational perturbation 
is reached when H becomes unity. Basically, the physical 
reason is the following: In vacuum the electromagnetic 
and gravitational modes obey the same dispersion rela- 
tion, and therefore - due to the mode coupling provided 
by the unperturbed magnetic field - the system evolves 
in a non-stationary way. In particular gravitational wave 
energy may be continuously converted into electromag- 
netic wave energy, as will be examined in more detail be- 
low. In the presence of a plasma, however, the induced 



2 



currents change the dispersion relation of the electromag- 
netic wave, and the resulting detuning of the modes sat- 
urate the conversion of energy between them, making a 
steady state solution (in a frame moving with the velocity 
of light) possible in principle. For a strongly magnetized 
plasma, on the other hand, the induced plasma currents 
cannot grow continuously with h, as we have seen above. 
For sufficiently high gravitational amplitude this means 
that the plasma currents are of little significance, prac- 
tically the plasma appears as vacuum for Ti. > 1, and in 
particular solutions depending only on z — t are impossi- 
ble. This conclusion is not dependent on the absence of 
thermal effects in our calculations in section II B. Gen- 
erally the addition of thermal motion only modifies our 
expressions ( |ll|) by a factor of the order l + {vt/c)^, where 
vt is th e thermal velocity. In particular, the divergence 
of ( lib ) still occurs for a finite value of H. 

On the other hand, it is clear that our omission of the 
back reaction of the electromagnetic wave on the gravita- 
tional pulse in principle could change this picture, since 
obviously certain components of the energy momentum 
tensor also diverges when H ^ 1, implying that the grav- 
itational wave amplitude could indeed be diminished due 
to the infiuence of the generated EM- wave. The effects of 
the selfconsistent gravitational field caused by the plasma 
perturbations are discussed in the Appendix, but will be 
omitted here since it turns out that the backreaction on 
the gravitational wave is negligible in the application to 
be discussed in this article. 

Since it is clear that for H > 1 the generated currents 
cannot stop the growth of the EM-wave, we simplify the 
picture from now on by putting the density to zero and 
thus totally ignoring the plasma effects. The general so- 
lution to Eq. (0) for the electric field in the presence of 
a monochromatic gravitational wave h = hcos[k{z — t)] 
can then be written 



6B = Ej_ = \k{C,,z + Ctt)Bohsin[k{z - t)] 

+E+{z -t) + E^{z -t) , 



(13) 



where Cz + Ct = I and E+ and E- are arbitrary func- 
tions. For an initial value problem where the plasma is 
unperturbed in the absence of the pulse Cz — 0, Ct — 1 
and i?+ — E- = 0, i.e. the electromagnetic amplitude 
grows linearly with time. For a boundary value problem, 
on the other hand, where the external magnetic field Bq 
occupies a region z > and there is a gravitational wave 
but no EM-waves propagating into the magnetized re- 
gion, clearly Cz = I, Ct = and E+ = E- = 0, i.e. we 
have a linear spatial growth instead. For the applications 
to be discussed later on we will be interested in a situa- 
tion where Bq is not necessarily static. We thus note that 
qualitatively the solution given by Eq. (13) still appHes 
for a quasi-static situation, i.e. where the dependence of 
Bq on time is slow enough such that the electric fields E 
associated with the time variations fulfills E/Bq ^ 1. 

In principle we can achieve very large EM-wave am- 
plitudes also when we abandon the specific solutions de- 



pending on z — t. However, since the growth is only linear 
in t and/or z, apparently we need large times/distances 
of coherent interaction. For a boundary value problem 
we can roughly define the effective distance of interaction 

as Zoff 



SBn 



ZcSBa^chaih'^ha 



(14) 



where the index char denotes the characteristic values 
of the various quantities in the region of interest, and 
the prime denote differentiation with respect to the ar- 
gument. 

To summarize: Eq. (lib) have a class of physically 



sound solutions, but also unphysical ones with the prop- 
erty Sn ^ oo as Ti ^ 1. The singular behavior is caused 
by the insistence to look for solutions that moves with a 
specific velocity, together with the omission of the self- 
consistent gravitational field from the plasma perturba- 
tions. The divergent solutions can be removed either 
by considering a boundary or an initial value problem, 
as discussed in this subsection, or by considering the 
backreaction of the plasma perturbations on the grav- 
itational wave, as discussed in the Appendix. The al- 
ternative considered here is the most relevant one with 
regard to astrophysical appHcations. Real astrophysical 
systems have finite distances of interaction between grav- 
itational waves and plasma waves, that can be estimated 
on physical grounds. Thus when estimating the maxi- 
mum magnetic field perturbation that can be produced 
by a gravitational wave i n a given situation, we can in 
principle apply solutions (lib) together with ( |T^ ) but we 
must note the upper bound for ^Bmax that exists for a 
given Zcff and is given by Eq. (F3). 



III. PHOTON FREQUENCY SHIFT 

We now consider the effect of the gravitational wave 
perturbations on high frequency photons in a plasma. 
For simplicity we assume that the photons propagate 
parallel to the gravitational waves and let them be rep- 
resented by the vector potential A — Aexp(i6') -I- c.c, 
where c.c. stands for complex conjugate. Making the 
approach of geometrical optics ||lO||, the wave number 
k = dzO and frequency lu = —dtO satisfies some local 
dispersion relation cu = W{z, t, k). The amplitude of the 
vector potential is assumed small and by high frequency 
photons we mean uj 3> Wp(i),ti;c(i)- 

Due to the gravitational waves the plasma has a back- 
ground of possibly large fields 5n, v\\, 6B, and E± all 
being functions of z — t and varying on a time and length 
scale much longer than that of A. 

Since uj ^ (^c{i) the high frequency pulse approxi- 
mately behaves as if the plasma is unmagnetized. The 
equation of motion linearized in the high frequency (hf) 
variables reads 



d 



d 



hf 



-g(») 



dA 
'dt 



dA' 



(15) 
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and thus v 



hi 



-(7(i)A/7m(j), where the large scale vari- 
ations have been neglected. The induced high frequency 
current is therefore j^^ = —ujpA/^Q, where the plasma 
frequency is Wp = (X^i '?fi)'^/^oW(i)7)^/^. Taking the time 
derivative of Ampere's law gives the following wave equa- 
tion for the photons 



9*2 







(16) 



We recognize the dispersion relation as oj = [k'^ +u!p{z — 
t)Y^^, where we assume that the variations in the plasma 
frequency are determined from Eqs. ( pj ) together with 

The change in the wave number and frequency as the 
wave propagates through the nonuniform and time vary- 
ing media with velocity Vg — duj/dk is given by the ray 
equations 



dk 
'dt 



dW 
dz 



'dt 



dW 



(17) 



Note that W is a function of z — * and introduce coor- 



dinates, ^ 



Vgt, T 



t locally moving with the pho- 



tons, i.e. it should be understood that Vn 



Vg{T = To) 



for some tq. Then, in a small neighborhood of tq it holds 
that du)/dT = -dW/d£,. Using 9^ = (1 - Vg)-^dr, this 
can be integrated from time 1 to 2 (which need not be a 
small interval), noting that 1 — « lJ^I2uP' . The result 
is 



CJ2 



(18) 



where the indices 1 and 2 denote the values at t\ and T2 , 
respectively. An interesting aspect of Eq. ([l^ ) is that the 
frequency conversion factor N — ujx /uj2 is independent of 
the frequency regime of the EM- wave. Thus, in princi- 
ple, x-rays can be turned into gamma rays, just as well as 
infra-red waves can be converted into the visible regime. 
This is in contrast to laser excited wake fields [H, where 
efficient frequency shifts can only take place provided the 
frequency of the converted pulse roughly lies in the same 
frequency regime as the exciting laser pulse. The reason 
for the difference is that the density gradients propagates 
with exactly the speed of light in our case, whereas, nat- 
urally, the corresponding velocity is slightly less than c 
in the laboratory experiments. The necessary distance 
of acceleration for a given conversion factor N is propor- 
tional to a;2, however, and this puts certain limits for the 
applicability to the highest frequency regimes. 



IV. EXAMPLE 

We have found that large density perturbations trav- 
ehng with the velocity of light can be induced by small 
gravitational wave perturbations, provided the cyclotron 



frequency is much larger than the plasma frequency, as 
described by Eqs. (|ll|). Furthermore, photons propagat- 
ing in a moving density gradient can undergo frequency 
up-conversion (or down-conversion), as described by Eq. 
(P^. In principle the effects can be large, even for a mod- 
erate deviation from fiat space-time. It is not yet clear 
that the predicted frequency conversion can be observed 
during reasonable conditions, however, and our aim in 
this section is to provide estimates to shed light on this 
question. In this section we reinstate the speed of light 
in all expressions. 

As a source of gravitational radiation we consider a bi- 
nary system. At least one of the objects should have a 
moderate to strong magnetic field (in order to make the 



parameter 



p(«) 



/^c{t)) small), and the objects should 



be compact (as to make the gravitational wave frequency 
and amplitude before merging large). Thus, for definite- 
ness (and calculational simplicity due to symmetries) we 
assume that the system consists of two neutron stars of 
equal mass M0, separated by a distance of 40 where 
Rs — 2GMq/c^ « 3 km. Furthermore, the surface mag- 
netic field of each neutron star is assumed to be 10® T. 
For the unperturbed plasma density profile, see Fig. 0. 

The area surrounding the binary system can loosely 
be divided into three regions (Fig. |l]). The interval 
20 Rs — 30 Rs from the center of mass (CM) roughly 
constitutes region I, which is the region where most of the 
gravitational energy is gained by the EM-wave. Using a 
Newtonian approximation, with d = aRs and r — (3Rs, 
it is straightforward to show that ~ (2q;/3)~^, where 
d is the separation distance between the binary objects 
and r is the observation distance from the center of mass 
of the system. In order to obtain an estimate of the am- 
plitude of the generated EM-wave, we combine the above 
expression for the gravitational wave amplitude with Eq. 
(P^ and the data given above. The result is 



5B_ 



7 X IQ-s 



(19) 



at the of end region I. In region II (approximate inter- 
val 30 Rs - 3500 Rs from the CM) 5B/Bo is still small, 
and - as seen by Eq. ( p^ ) - the relative density per- 
turbation is thereby small as well, which limits the fre- 
quency conversion effect in this region. However, the 
gravitationally induced EM-wave suffers spherical atten- 
uation, whereas the unperturbed magnetic field is that 
of a dipole, and consequently the relative density per- 
turbation grows quadratically with distance. The end 
of region II is defined as the necessary distance to make 
SB/Bq of the order unity due to this increase. (For 
pulsars with period longer than 35 ms, region I and II 
lies in the near zone, and thus the unperturbed magnetic 
field indeed decays cubically in the region of interest, al- 
though the unperturbed field becomes a radiation field 
outside the light cylinder of the pulsar.) In region III 
(approximate interval 3500 Rs — 10® Rs), the relative 
density perturbation is appreciable, and thus the main 
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frequency conversion occurs here [fT^ . 

At the beginning of region III the relative density per- 
turbation is Sn/no ~ 1, in agreement with Eq. ( F2| ). An 
EM-wave with initial frequency uj = Wmin = 10 rad/s 
can move from a density minimum to a density maximum 
during a "laboratory system distance" Lfrcq = cTficq = 
cigrad/(c - Vg) ~ w^ax-^grad/t^p, whcrc Lgrad is a typ- 
ical density gradient scale length. For definiteness we 
assume that the pulsars have periods of the order of 350 
ms, in which case SB/Bq may increase to 5B/Bq ~ 10 
for the most of region III. In our example the maximum 
frequency magnification N thus is 



N 



p,max 



p.min 



10 . 



(20) 



Inserting 



10^ rad/s, and letting 



p,max 



10 rad/s (corresponding to no ~ 10 cm~^) we obtain 
Lfrcq — 10^i?S; i-e. the acceleration can take place within 
region III. Strictly applying our one-dimensional calcula- 
tions of Sec. 01 means that frequency up-converted EM- 
waves will be down-converted and vice versa, if the gravi- 
tational source and the induced density perturbation are 
indeed periodic. In our example, on the other hand, the 
successive frequency conversion effects will decrease with 
the distance from the source, and thus for an earth based 
observer the radiation generated in region III should show 
periodic up- and down-conversions. The frequency con- 
version ratio of Eq. ( pO| ) is of course a maximum value 
of our example, that occurs for radiation generated at a 
density extremum, but all radiation generated in region 
III will be up- or down converted with a factor in the 
interval 1 — N, and consequently the effect should be ob- 
servable provided the object is close enough for radiation 
generated in region III, in the approximate frequency in- 
terval 10^^ rad/s < lu < 10^"* rad/s, to be detected, where 
the upper limit is imposed by the fact that the system has 
a finite distance of interaction. If we try to increase the 
interaction efficiency by considering higher plasma densi- 
ties the electromagnetic wave damping due to Thomson 
scattering becomes prohibitive ||l2|. 



V. SUMMARY AND DISCUSSION 

We have considered the generation of traveling den- 
sity perturbations in a magnetized plasma induced by 
gravitational radiation. Provided X]i(^p(i)/'^c(i)) ^ 
significant density perturbations, i.e. Sn/no ~ 1, can be 
induced even by a small gravitational wave with h <^ 1, 
provided 7i ~ 1. Basically the large effect is possible be- 
cause of the approximate agreement of the dispersion re- 
lations between the fast magnetosonic and gravitational 
modes in the regime J2ii^p(i)/^c{i)) ^ ^' which in turn 
allows for a long distance of coherent interaction. 

In order to find a mechanism where the induced density 
perturbations may give rise to earth-based observational 



effects, we have studied frequency conversion of electro- 
magnetic wave packets traveling in the moving density 
gradients. The formula (18), relating the frequency of 
the wave packet for two different positions in the moving 
density profile, is in conceptual agreement with the cor- 
responding results of Ref. |^ , which considered an anal- 
ogous situation but where the density perturbation was 
due to plasma oscillations traveling with a phase veloc- 
ity slightly less than the speed of light c. In our case 
the gradients move with exactly c, however, and thereby 
the maximum frequency conversion factor N does not de- 
crease with the initial frequency [as for conventional pho- 
ton acceleration], in principle allowing for up-conversion 
even of 7-rays. 

The ideaHzations made in Sees. II and |l| is a some- 
what too strong for our results to be directly applicable to 
a situation of astrophysical relevance. In particular, we 
cannot consider the unperturbed plasma as homogeneous 
and the geometry as one-dimensional when making esti- 
mates. In our example with a binary system as a source 
of gravitational radiation, we have thus been forced to di- 
vide the neighborhood of the system into three regions: 
Region I where most of the energy transfer into electro- 
magnetic wave energy occurs, region II where the relative 
density perturbation grows, and region III where the fre- 
quency conversion takes place. In order to describe the 
physics in region I adequately we must abandon solu- 
tions that depend on z — ct only , and the basis for this 



has been discussed in Sec. [IC. By making estimates 



based on our analytical calculations, we conclude that 
the gravitational waves emitted by a system of binary 
pulsars close to merging may result in periodic frequency 
up- and down-conversions of electromagnetic radiation 
in the infrared part of the spectrum. The frequency of 
the up- and down-conversions coincides with the gravita- 
tional wave frequency, i.e. it is twice the orbital frequency. 



VI. APPENDIX 

In this appendix we are going to investigate the regime 
of validity for the multi-component test fluid approach. 
Normally we think that by continually decreasing the pa- 
rameters proportional to the unperturbed energy density, 
at some point the fluid in an external gravitational fleld 
can be treated as a test fluid. In our case the situation 
is not quite that simple, since we can decrease the elec- 
tromagnetic (oc Bq) and the rest mass energy density 

'l{^)/'^c{^)) constant. 



(cx no)at the same rate keeping 
Since our solution in section II B has a diverging energy- 
momentum tensor whenever Ti = '^'h/ J2ii'^p(i)/'^c{i)) ~^ 
1, clearly we cannot justify the test fluid approach sim- 
ply by assuming a sufficiently low unperturbed energy 
density. To shed light on the physical effects due to the 
selfconsistent gravitational field, we will first consider the 
linearized theory. This will provide a guide for making 
estimates of the regime of validity of our (nonlinear) test 
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matter solution in section II B, and also makes it possi- 
ble to justify the omission of selfconsistent gravitational 
effects in section IV . 

We divide all quantities into an unperturbed part (i.e. 
the value in the absence of the gravitational perturba- 
tion) and a perturbed part. We note that the only vari- 
ables that are nonzero in the unperturbed state are the 
density (= no), the magnetic field {— -BqCi) and the met- 
ric (= Tjf^^). It should be emphasized that in addition to 
the direct effect on the dispersion relation from the mat- 
ter, which we will study below, there is also an indirect 
contribution (that will be omitted here) to the dispersion 
relation from the background curvature produced by the 
(unperturbed) matter. In the regime where the gravita- 
tional wave length is much shorter than the background 
curvature, however, the shortwave approximation can be 
applied, which imply that these two effects can be stud- 
ied separately and their contribution to the dispersion 
relation of the gravitational wave can be added, see e.g. 
Ref. . In the above scenario (provided thermal effects 
are still neglected) the only effects from the gravitational 
wave on the plasma perturbations are from the effective 
currents in (|4a|)-(4b), where, in the present case, we have 

= —j^ = —{l/2)Bo{dh/dz) and the other components 
are zero. Thus using Maxwells eqs. and the the set of 
fluid equations for each particle species and the same ap- 
proximations as in section II (but avoiding the ansatz 
d/dt — d/dz) we will obtain a wave equation for the fast 
magnetosonic wave, modifled from the standard textbook 
form by allowing for an arbitrary value oiYliii^p^i)/ 
and with a gravitational "source term" due to the effective 
gravitational currents above. The result is 



r2 



92 



1 + C\ 



5B 



-^Bo (21) 



where we have introduced the Alfven velocity C'a — 

( Z]i(^p(i)/^c(i))) "(Note that Ca may be larger than 
unity, but, as can be seen above, the actual magnetosonic 
wave velocity is smaller or equal to C^-) The system 
is closed selfconsistently by Einstein's fleld equations, 
which, after linearization reduces to (cf Eq. 4.9 in Ref. 

@) 



92 



722]lin — — BqSB 

Mo 



(22) 

where lin stands for "linear part of". It is simple to com- 
bine Eqs. (^ and ( p2| ) into a single wave equation for 
the coupled fast magnetosonic and gravitational mode. 
However, it is probably more illustrative to proceed by 
considering the corresponding dispersion relation. Mak- 
ing a plane wave ansatz, SB = (5i?exp[i(fcz — ujt)] and 
h — hexp[i{kz ~ ujt)], we directly flnd the dispersion re- 
lation: 



Z2itGBI 
Mo 



.2-PCl/(l + Cl) 



(23) 



from Eqs. (|l|) and (p2|). Thus the presence of matter 
causes a phase velocity u/k > 1 and a group velocity 
dcu/dk < 1. A further consequence is that the gravi- 
tational wave also becomes dispersive. Apparently the 
relation between SB and h is 



SB = Bah 



k^cy{i + ci) 



(24) 



where the omission of the selfconsistent gravitational field 
is a valid approximation only if we can use the vacuum 
dispersion relation uj'^ — k'^ = as an approximation in- 
stead of Eq. ( p3| ) when calculating SB from ([24|). From 
now on we will focus on the regime ^ 1 , which make 
the magnetosonic phase velocity close to unity. Since the 
(typically small) right hand side of ( |2^ ) now must be com- 
pared to the small phase velocity difference of the (un- 
coupled) magnetosonic and gravitational waves, the con- 
dition for omitting the selfconsistent gravitational field is 
significantly stronger if one should get an approximately 
correct magnetic field, and not just a small contribution 
from the right hand side in the dispersion relation (p3|). 
For Ca ^ 1 the condition for omitting the selfconsistent 
gravitational field, and still obtaining an approximate ex- 
pression for SB, becomes: 





Mo C\ 



(25) 



The above validity condition is obtained by comparing 
the magnetic field obtained from the full selfconsistent 
dispersion relation and its vacuum approximation. A 
much simpler way to arrive at the same condition as in 
( p5| ) is to demand that the relative contribution from the 
energy momentum tensor terms in Einstein's equations 
should be much smaller than the relative velocity differ- 
ence between the magnetosonic and gravitational waves. 
The advantage with this latter formulation of the valid- 
ity condition is that it can be easily applied also when 
the relation between SB and h as well as the expression 
for the energy momentum tensor are nonlinear. Adopt- 
ing this condition for omitting the selfconsistent grav- 
itational field when the plasma response to the metric 
perturbation is nonlinear we write 



327rGmax((5r) < 



char : 
^A 



^cha 



(26) 



where max{ST) denotes the maximum deviation from the 
unperturbed value of the perturbed energy momentum 
tensor for any of its components, and the index "char" 
denotes the characteristic value of the gravitational wave 
frequency and metric perturbation, respectively. For the 
regime when Eq. ( |2^ ) is violated, obviously our solution 
in section II B must be modified to take the selfconsistent 
gravitational field into account, and this may result in 
new types of solutions describing, for example, nonlinear 
solitary gravitational pulses. This problem is outside the 
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scope of our present article, however. We note that our 
example in section IV, fulfills the validity condition ( |2^ ) 
with a margin of several orders of magnitude. 
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FIG. 1. The neighborhood of the binary system is di- 
vided into three regions: Region I (20iis — 30iis), region 
II (SOiis - 3500J?s), and region III (3500J?s - lO'' i?s). In 
regions I and II we are situated in the near zone of the mag- 
netic field of the pulsar. Thus the plasma density is low, and 
we assume that the plasma particles do not interfere with the 
approximations made in the example. For this to be true, 
7i S> 1 must hold in region I, which is satisfied even for very 
high densities. Furthermore X!i('^p(i)/'^c(i)) ^ 1 should apply 
in region II, which is fulfilled for no < 10^ cm~"^. In region III, 
which is mainly outside the light cylinder of the pulsar, we 
assume the plasma density no to be of the order of 10^^ cm~^. 
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